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Abstract 

Starting from the Hamiltonian formulation for the inhomogeneous Mixmaster dynam- 
ics, we approach its quantum features through the link of the quasi-classical limit. We fix 
the proper operator-ordering which ensures that the WKB continuity equation overlaps 
the Liouville theorem as restricted to the configuration space. We describe the full quan- 
tum dynamics of the model in some details, providing a characterization of the (discrete) 
spectrum with analytic expressions for the limit of high occupation number. One of the 
main achievements of our analysis relies on the description of the ground state morphol- 
ogy, showing how it is characterized by a non-vanishing zero-point energy associated to 
the Universe anisotropy degrees of freedom 

PACS: 04.20.Jb, 98.80.Dr, 83C 

1 Introduction 



Einstein's equations, despite their non-linearity, admit a generic solution asymptotically to the 
cosmological singularity having a piece-wise analytic expression. This dynamical regime was 
derived by V. A. Belinski, I. M. Khalatnikov and E. M. Lifshitz (BKL)pQ [2] and it has an 
oscillatory-like behavior which extends point by point in space the homogeneous Mixmaster 
dynamics of the Bianchi type IX model [3] . When a generic initial condition is evolved toward 
the initial singularity, a space-time foampfl, [5] appears: this is a direct consequence of the dy- 
namical decoupling which characterizes space points close enough to the singularity and which 
can be appropriately described by the Mixmaster point like measure HE] (a wide literature 
exists about the chaotic properties of the homogeneous Mixmaster model; see, for instance, 
[H El El HO] and the references therein) . 

However, this classical description is in conflict with the requirement of a quantum behavior 
of the Universe through the Planck era; in fact there are reliable indications [8] that the Mix- 
master dynamics overlaps the quantum Universe evolution, requiring an appropriate analysis 
of the transition between these two different regimes. 

Indeed the dynamics of the very early Universe corresponds to a very peculiar situation with 
respect to the link existing between the classical and quantum regimes because the expansion 
of the Universe is the crucial phenomenon which maps into each other these two stages of the 



evolution. As shown in [8] the appearance of a classical background takes place essentially 
at the end of the Mixmaster phase when the anisotropy degrees of freedom can be treated as 
small perturbations; this result indicates that the oscillatory regime takes place almost during 
the Planck era and therefore it is a problem of quantum dynamics. However the end of the 
Mixmaster (and in principle the quantum to classical transition phase) is fixed by the initial 
conditions on the system and in particular it takes place when the cosmological horizon reaches 
the inhomogeneity scale of the model; therefore the question of an appropriate treatment for 
the semiclassical behavior arises when the inhomogeneity scale is so larger than the Planck scale 
that the horizon can approach it only in the classical limit. Some interesting features of the 
quantum Mixmaster have been developed by [TT], [12], [13l [14] and [T5l [16] for the homogeneous 
and the inhomogeneous cases, respectively. 

Here we refine this analysis both connecting some properties of the quantum behavior to the 
ensemble representation of the model, both describing the precise effect of the boundary con- 
ditions on the structure of the quantum states. 

We start by an Arnowitt-Deser-Misner (ADM) [TT] representation of the system dynamics which 
allows us to disregard the contributions of the spatial gradients relative to the configurational 
variables , thus reducing the dynamics to a number of oo 3 independent point-like Mixmaster 
model. Such representation can take place independently of the gauge choice [18], nevertheless 
here we require one of the Misner-Chitre like variables to play the role of time for the system. 
This way, we can describe the whole evolution in terms of a triangular billiard on the Poincare 
upper half-plane where the point Universe randomizes. 

Then we apply the Liouville theorem restricted to the configuration space using the Hamilton- 
Jacobi solution of the model, and and we require the continuity equation to match the WKB 
limit. Such a requirement leads us to determine the proper operator-ordering associated to the 
ADM operator in a unique way. We emphasize that the existence of an energy-like constant of 
motion not only provides a microcanonical measure for the statistical dynamics but also induces 
a quantum dynamics completely described by a time-independent Schrodinger equation on the 
Poincare plane. However, the effective Hamiltonian associated to such an eigenvalue problem 
is non-local, characterized by the presence of a square root function; in agreement with the 
analysis developed by [19], we make the well-grounded hypothesis that the eigenfunctions form 
is independent of the presence of the square root, since its removal implies the square of the 
eigenvalues only. 

A crucial step relies on recognizing how our squared equation can be restated as the Laplace- 
Beltrami eigenvalue problem. Hence a characterization of the "energy" spectrum comes out 
from imposing the triangular boundary conditions (the semicircle is replaced by a straight line 
preserving the measure of the domain). 

We outline that the spectrum is discrete and that admits a point-zero energy which implies an 
intrinsic anisotropy of the quantum Mixmaster Universe. We numerically investigate the first 
"energy" levels and finally provide an analytic expressions for high occupation numbers. 
The structure of the presentation is the following: Section II is devoted to fix the classical 
Hamiltonian formulation of the system dynamics, based on the use of Misner-Chitre-like vari- 
ables which allow the ADM reduction of the degrees of freedom. In Section III we fix the 
solution of the Hamilton- Jacobi equation within the domain allowed by the vanishing of the 
potential. Such a solution is at the ground level of the analysis presented in Section IV; in fact 
here we restrict the Liouville theorem to the configuration space, eliminating the momentum 
dependence by virtue of an integration along the Hamilton- Jacobi dynamics. In Section V we 
describe the Schrodinger equation for the Mixmaster model and construct the WKB limit; by 
the comparison of this semiclassical behavior with the equation induced by the Liouville theo- 
rem, we get a unique choice for the operator-ordering of the super-Hamiltonian kinetic term. 
Finally, in Section VI we provide a detailed description of the Mixmaster spectrum and eigen- 
functions; particular attention is devoted to the ground state in view of the idea that it is the 



expected Universe configuration during the Planck era. Section VII is devoted to concluding 
remarks on the physical issues of our discussion, and particular attention is devoted to fix the 
full inhomogeneous quantum picture. 

Over the whole paper we adopt units such that c = 16nG = 1. 



2 Hamiltonian Formulation 

In this section we introduce the Hamiltonian formalism for a generic inhomogeneous model in 
vacuum in terms of the Misner-Chitre-like variables [3| 120]. 

In the ADM approach, the line element associated to such a model can be written in the form: 

ds 2 = N 2 drf - -f a p{dx a + N a drj)(dx p + JV%) (1) 

where N and N a denote the lapse function and the shift- vector, respectively, while 7 a/ 3 (a, (3 = 
1, 2, 3) is the 3-metric tensor of the spatial hyper-surfaces S 3 n = const, and we take [4] 

1ai3 = e qa 5 ad O a b O d c ^^, a, b, c,d,a,P = 1,2, 3, (2) 

where q a = q a (x,n) and y b = y b (x,n) are six scalar functions, and O b = O b (x) a SO (3) matrix 
(repeated indexes are summed). 

Using the metric tensor ([2]), the action for the gravitational field reads 

S = [ d7]d 3 x (p a d v q a + U d d v y d - NH - N a H a ) , (3) 



with 

h 1 

'1 



(4) 



H a = U c d a y c + p a d a q a + 2 Pa (0- 1 ) b a d a O a b ; (5) 

H and H a being the super-Hamiltonian and the super-momentum, respectively. 

In 03]) and (jSJ) p a and Il c play the role of conjugate momenta to the variables q a and y b , 

respectively, and is the Ricci 3-scalar which behaves like a potential term. 

No sooner should we adopt y a like new spatial variables, than the Super-momentum constraint 

H a — (provided by the variation of ([3]) with respect to N a ) can be solved, and the action 

rewrites as[T%] 

S = [ d V d 3 y (p a d,q a + 2p a (0~ 1 ) c a d v O a c - NH) . (6) 

Furthermore the potential term appearing in the Super-Hamiltonian (J3J) can be approximated 
toward the singularity as 

U = J2®(Qa)i (7) 

a 

where 

This picture arises from the vanishing of the metric tensor determinant close to the singularity, 

fmH thp mifmtitips O _ s known in thp litprntnrp fic: nniftnt.rnrtu nnrnrnpfprs put a pln^prl Hnmain 



Tq to which the dynamics is restricted. 

By virtue of the potential structure (ED), the Universe dynamics evolves independently in each 
space point inducing a corresponding factorization for the phase space of the model. 
Since in Tq the potential U asymptotically vanishes, close to the singularity the relation 
dpa/drj = holds and then the term 2p a {0~ 1 ) c a d r) 0'^ in ([6]) behaves like an exact time-derivative 
and can be ruled out of the variational principle. 

The ADM reduction is completed by introducing the so-called Misner-Chitre-like variables 
[H HQl 1201 |2T] .whose relevance consists in making the anisotropy parameters independent of r, 
which will behave as the time variable: 



q 



q = e 
q 3 



V^ 2 - l(cos 9 + v/3 sin 9) - £ 
V^ 2 - l(cos0- v^sintf) - i 
-e T + 2v/f 2 - lcosfl 



In terms of these new variables the super-Hamiltonian rewrites 



(10) 



(11) 



and the Q a 's become 



Qi 
Q2 
Q3 



1 



1 



3 3£ 

1 



cos 



9 + ^3 si 



sm( 



3 3e ; 

1 , 2^e^i 



+ 



cos 
cos 9 



(12) 



Let's solve the constraint H = (obtained variating N in Q) with respect to p T in the domain 
Tq in order to perform the ADM reduction, thus obtaining: 



-Pr = e 



(e - i)pI + 



(13) 



Then, taking the time gauge d v r = 1, the reduced action explicitly reads as 



drd 3 y (p^d T ( + pgd T 9 



(14) 



As we approach the singularity, e behaves like a constant of motion, i.e. de/dr = de/dr = =>- 
e = E{y a ). 

The dynamics of such a model is equivalent to (the one of) a billiard-ball on a Lobatchevsky 
plane [Tl [Til [T5| l2"2"]; this can be shown by the use of the Jacobi metric, i.e. a geometric approach 
that reduces the dynamical equations of a generic system to a geodesic problem on a manifold. 
Such a technique applied to (fTI]) produces a geodesic equation corresponding to the line element 



di 1 = E 2 {y a ) 



de 



(e - 1) 



(15) 



The manifold described by ([TBI) turns out to have a constant negative curvature, where the 
Ricci scalar is given by R = —2/E 2 : the complex dynamics of the generic inhomogeneous 

mnHpl rpcniltc: in a nollpptinn nf HppnnnlpH Hvrmmipfil swtpms nnp for Pfiph nnint nf thp snapp 
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Figure 1: The domain Tq(u,v) in the Poincare upper half-plane where the dynamics of the 
space-point Universe is asymptotically restricted by the potential. 
The measure is finite and equal to it. 



and all of them equivalent to a billiard problem on a Lobatchevsky plane. 

Among the possible representations for it, we choose the so-called Poincare model in the complex 

upper half-plane [6] that can be introduced with the following coordinate transformation 

_ 1 + u + u 2 + v 2 



^ V r, , (16) 



9 = — arctan - 

-1 + 2u + 2u 2 + 2v 2 

The line element for this 2-dimensional surface now reads 

, 2 du 2 + dv 2 . . 

ds 2 = - 17 

v 

Figure (pQ) shows how, in terms of these coordinates, the anisotropy functions cut a stripe in 
the plane (u, 0, v) characterized by a finite measure /i = 7r and three unstable directions to 
asymptotically escape. 

The boundary is given by the set of points where one of the anisotropy parameters is equal to 
zero; it is a geodesic triangle, where the edges are given by the vanishing of the numerators of 
QaS, now reading 

Qi(u,v) = —u/8 

Q 2 (u,v) = (l + u)/5 (18) 
Q 3 (u,v) = (u 2 + u + v 2 )/5 

S = u 2 + u + l + v 2 
In the (u, v) scheme, the ADM Hamiltonian ffl3l) assumes the expression 

Hadm = e = v^/pl+pl, (19) 
which will be the starting point of our analysis. 

We conclude by observing how [TJ [2] the present representation of the Lobatchevsky plane 
is suitable to link this Hamiltonian formulation to the piece-wise representation provided by 
Belinski, Khalatnikov and Lifshitz map, as it is easy to recognize that for v = 0, the functions 
Q a reduce to the familiar Kasner indexes. 



3 Hamilton-Jacobi Approach 



We shall now derive the Hamilton-Jacobi (HJ) [23] equations for the system (fl9|) in view of the 
following developments. 

The HJ prescribes to change the momentum variables into the derivatives of a functional S 
with respect to the configurational variable: implementing the HJ technique, the Hamiltonian 
relation ffl9l) leads to the functional differential equation in Tq 



« J( s 4)\( s 4) 2 . (20) 



5t y \5u J \ Sv J 
To obtain the solution of this dynamical equation, we take S in the form 



S = S (u,v)-[ J d A ye{y a ) ) r, (21) 

where So satisfies the equation 
and it is therefore provided by 



5S \ 2 , / 5S ^ " 



S (u, v) = k(y a )u + v/ e 2 _ k \y")v* - e In ( 2 t + VfZE^l \ + c (^) ; ( 23 ) 



e 2 v 



where k(y a ),c(y a ) are integration constants. 

The expression (BIT) , together with (1231) and the features of the potential wall (|T8|) . summarizes 
the classical dynamics of a generic inhomogeneous Universe. 



4 The Statistical Mechanic Description 

In the previous Section 2 we have pointed out that, for the Mixmaster inhomogeneous dynamics, 
the spatial points decouple approaching the singularity and an energy-like constant of motion 
appears; let us discuss the problem from a statistical mechanics point of view by treating the 
system microcanonical ensemble. 

The physical properties of a stationary ensemble are described by a distribution function p = 
p(u,v,p u ,p v ), obeying the continuity equation defined in the phase space (u,v,p u ,p v ) as 



du ^ dv 
where the dot denotes the time derivative 
du dH ADM v 2 



d(up) d{vp) d{p u p) d(p v p) 



+ 



dp u 



+ 



dp v 







(24) 



u 



r = 



dr 
dv 



dp u 
dH ADM 

dp v 



-Pu 



-Pv 



Pu = 

Pv = 



dp u 
dr 
dp v 
dr 



dH 



ADM 



du 
dH ADM 



dv 



e 
v 



(25) 



We stress how the above continuity equation provides an appropriate representation for the 
ensemble associated to the Mixmaster only when we are sufficiently close to the initial singu- 
larity and therefore the infinite-potential wall approximation works; in fact, such a model for 

thp nntpntifil tprm nnrrpsnnnrls tn HpfiI with thp pnprcrv-likp pnnRtfint nf thp motion Finrl Iiyps 



the microcanonical nature of the ensemble. From a dynamical point of view this picture arises 
naturally because the Universe volume element vanishes monotonically (for non stationary cor- 
rection to this scheme in the Misner-Chitre like variables see [24|). 

We are interested in studying the distribution function in the (u, v) space, and thus we will 
reduce the dependence on the momentum variables by integrating p(u, v, p u , p v ) in the momenta 
space and using the informations contained in the HJ approach. If we assume p to be a regular, 
vanishing at the infinity of the phase space, limited function of its arguments, and, if use the 
([231 . we work out the following equation for w(u, v; k) = f p{u, v , p u , p v )dp u dp v 



dw(u,v;k) f E\ 2 ^dw(u,v;k) E 2 — 2k 2 v 2 w(u,v;k) . . 

du y V W dv + kv 2 y/E 2 - (kvf ~ 

admitting a solution in terms of a generic function g of the form 



9[u + Vy/-^x-l 

w(u,v) = ^ — 7 - (27) 

vyE 2 — k 2 v 2 

The distribution function in (u, v) is obtained after the integration over the constant k. Indeed, 
this constant expresses the freedom of choosing the initial conditions (125]) . which cannot affect 
the chaotic properties of the model. Therefore we define the reduced distribution w(u,v) as 

w(u, v) = / w(u,v;k)dk, (28) 
J A 



where the integration is over the classical available domain for p u = k 



A 



E E 

v ' v 



(29) 



It is easy to verify that the microcanonical Liouville measure [H EJ [7J [10l [HI [2H [25] w mc (after 
integration over the admissible values of e ) corresponds to the case g = const, i.e. we get the 
normalized distribution 

E_ 

/ v 1 7T 
dk = — . (30) 

Summarizing, we have derived the generic expression of the distribution function for our model, 
fixing its form for the microcanonical ensemble which, in view of the energy-like constant of 
the motion e, is the most appropriate to describe the Mixmaster system when the picture is 
restricted to the configuration space. 



5 Quasi-classical Limit of the Quantum Regime 

Let us underline some common features between the classical and the semi-classical dynamics 
with the aim of fixing the proper operator-ordering [26] in treating the quantum approach. 
In fact, considering the WKB limit for h — > 0, the coincidence between the classical distribution 
function w mc (u,v) and the quasi-classical probability function r(u, v) takes place for a precise 
choice of the operator-ordering [16] only (for a connected topic see [8] and for a general discussion 
see [27]). 

As soon as we implement the canonical variables into operators, as in the canonical quantization, 



i.e 



Pa 



-ih- 



Pr 



-ih- 



OV OU OT 

the quantum dynamics for the state function <&(u,v,t) obeys, in each point of the space, the 
Scrhodinger equation 



ih 



d<&(u, v, t) 



H adm $(u,v,t) = h\ 



du 2 



,2-0 



dv 



dv J 



$(tt, v, t) 



(31) 



where we have adopted a generic operator-ordering for the position and momentum operators 
parametrized by the constant a |28j . 

In the above equation the Hamiltonian operator has a non-local character as a consequence 
of the square root function; in principle this constitutes a subtle question about the quantum 
implementation though. We will make the well-grounded ansatz[19] that the operators Hadm 
and H\ DM have the same set of eigenf unctions with eigenvalues E and E 2 , respective!;^. 



If we take the following integral representation for the wave function $ 

/oo 
^(u,v 1 E)e- lET/h dE, 
-oo 

the eigenvalues problem reduces to 



(32) 



H 2 ^(u,v,E) = ft 2 



— v 



2-a 



d_ 

dv 



dv 



V(u,v,E) = E 2 V(u,v,E) 



(33) 



In order to study the WKB limit of equation (|33l) . we separate the wave function into its phase 
and modulus 

*(u, v, E) = y/r(u,v,E)e itT( - u > v > E)/h . (34) 

In this scheme the function r(u,v) represents the probability density, and the quasi-classical 
regime appears as we take the limit for h — > 0; substituting ( 1341) in ( 1331) and retaining only the 
lower order in h, we obtain the system 



fdaV + fda\ 
\du) \ dv J 
dr da dr da 



du du dv dv 



E 2 



a da d 2 a d' 



a 



v dv dv 2 du 2 



(35) 



0. 



In view of the HJ equation and of the Hamiltonian ( 1191) . we can identify the phase a to the 
functional So of the HJ approach. 

Now we turn our attention to the equation for r(u,v); taking (123]) into account it reduces to 



k 



dr(u, v, E) 
du 



u9 dr(u,v,E) a{E 2 - k 2 v 2 ) - E 2 
— k l 1 ; == — r{u, v, E) = . (36) 



dv 



WE 2 - k 2 v 2 



Comparing (136!) with (|26l) . we easily see that they coincide (as expected) for a = 2 only. 
It is worth noting that this correspondence is expectable once a suitable choice for the con- 
figurational variables is taken; however here it is remarkable that it arises only if the above 
operator-ordering is addressed. It is just in this result the importance of this correspondence 
whose request fixes a particular quantum dynamics for the system. 



x The problems discussed in this respect by [29 do not arise here because in the domain Tq our ADM 
Hamiltonian has a positive sign (the potential vanishes asymptotically) 



Summarizing, we have demonstrated from our study that it is possible to get a WKB corre- 
spondence between the quasi-classical regime and the ensemble dynamics in the configuration 
space, and we provided the operator-ordering when quantizing the inhomogeneous Mixmaster 
model to be [28J 

^ - -4 M) • < 37 > 

By means of these results, we now face the full quantum dynamical approach. 



6 Canonical Quantization and the Energy Spectrum 

Our starting point is the point-like eigenvalue equation (133]) (for a discussion of this same eigen- 
value problem in the Misner variables see [3]) which, together with the boundary conditions, 
completely describes the quantum features of the model, i. e. 

du 2 dv 2 dv \ h 



V(u,v,E) = 
V(dT Q ) = 0. (38) 



In equation ([381 we can recognize a well known operator: by redefining ^f(u, v, E) = if)(u, v, E)/v, 
we can reduce (138]) to the eigenvalue problem for the Laplace-Beltrami operator in the Poincare 
plane 

^ + ^)i(>(u,v,E) = E a 'tl>(u,v,E), (39) 

which is central in the harmonic analysis on symmetric spaces and has been widely investigated 
in terms of its invariance under SL(2, C)(for a detailed discussion, see [30] and the bibliography 
therein). 

Its eigenstates and eigenvalues are described as 

ip s (u,v) = av s + bv l ~ s + v /w^ n ^o a «^-i/2(27r|n|t;)e 2 ™, a, b, a n G C 

V lb ^j s (u,v) = s(s-l)ip s (u,v) (40) 

where K s _ 1 / 2 {'27rnv ) are the modified Bessel functions of the third kind, and s denotes the index 
of the eigenf unction. This is a continuous spectrum, and the summation runs over every real 
value of n. 

The eigenf unctions for our model read as 

v, E) = av*- 1 + bo'* + a n Ks - 1/2{ ^ lnlv) e^ nu (41) 

with eigenvalue 

E 2 = 8 (1 - s)h 2 . (42) 

The spectrum and the explicit eigenfunctions are obtained by imposing the boundary condi- 
tions (1351) . i.e. requiring that equation (T4T1) vanishes on the edges of the geodesic triangle of 
FigtH 

Since from our analysis no way arose to impose exact boundary conditions we approximate the 
domain with the simpler one in Figj2j the value for the horizontal line y = 1/n provides the 
same measure for the exact and the approximate domain. 
The conditions on the vertical lines u = 0, u = — 1 require to disregard the first two terms in 



Figure 2: The approximate domain where we impose the boundary conditions. 

The choice v — 1/tt for the straight line preserves the measure \i — 7T 



(1411) (a = b = 0); furthermore we get the condition on the last term 



CO 





(43) 



n>0 



which is satisfied by requiring ir s _!/ 2 (2n) = only. 

The functions K v {x) are real and positive for real argument and real index, therefore the index 
must be imaginary, i.e s — - + it. In this case these functions have (only) real zeros, and the 
corresponding eigenvalues turn out to be real and positive 



We remark that our eigenf unctions naturally vanish as infinite values of v are approached. 
The conditions ( 1431) cannot be solved analytically for all the values of n and t, and the roots 
must be worked out numerically for each n; there are several results on their distribution that 
allow us to find at least the first levels. A theorem on the zeros of these functions states that 
K iy {yx) = 0-v^0<x<l (for a proof see [33J); by this theorem and the monotonic dependence 
of the energy fT44|) on the zeros, we can easily search the lowest levels by solving equation (1431) 
for the first n. In Figj3]and Tab [1] we plot the first roots, and in FigJUand in Tabj2]we list the 



Below we will provide an analytical treatment of the high energy levels associated to our 
operator in correspondence to certain region of the plain (t, 0, n). In [13] the structure of the 
high energy levels is also connected to the so-called quantum chaos of the Mixmaster model, 
which was studied from the wave functional point of view in [12J and from the the path integral 
one in [T3] using Misner variables. Our analysis implicitly contains the information about the 
quantum chaos of the considered dynamics; in fact we can take a generic state of the system 
£(t, u, v) in the form 



{E/hf =t 2 + 1/4 . 



(44) 




£(t, u, v) = / dt c t , n ^n(r, u, V 



)5(K tt (2n))e- 1 ^ 



T 



(45) 



n 



2 for a detailed numerical investigation of the energy spectrum of the standard Laplace-Beltrami operator, 
especially with respect to the high-energy levels, see[3Tl[32], where it is also numerically analyzed the effects on 
the level spacing of deforming the circular boundary condition toward the straight line 



n 




10 



12 



Figure 3: The intersections between the straight lines and the curves represents the roots of 
the equation Kn{n) = 0, where K is the modified Bessel function 



n=2 


n=4 


ri=6 


n=8 


4.425 


7.016 


9.434 


11.768 


6.333 


9.353 


12.076 


14.655 


7.947 


11.313 


14.283 


17.059 


9.410 


16.264 


16.263 


19.212 


10.774 


17.742 


18.096 


21.203 



Table 1: We numerically investigated the roots of K it (n) = with respect t for different values 
of n. The values in the columns are the first t's in the above figure. 



and its evolution from a generic initial condition ^o(u, v) = £(to, u, v) at an initial instant tq 
provides all the quantum properties of the system. The quantum chaos is recognized in [12] by 
numerically integrating the Wheeler- De Witt equation from a gaussian like initial packet and 
outlining the appearance of a fractal structure in the profile of the resulting wave function; in 
our approach (as the infinite potential walls approximation works) the dynamics is provided by 
(145!) and evolving it from an initial localized wave packet the quantum chaos has to arise. 
About the information on the quantum chaos emerging from the high-energy spectrum we 
emphasize the following two important points: i) the non- stationary corrections due to the 
real potential term are expected to be simply small perturbation to our result rapidly decaying 
as the singularity is approached; ii) the analytic expressions we are going to provide for large 
values of t cannot be used to fix the existence of the quantum chaos because they explore 
limited regions of the plain (t, 0, n) only, but they are very useful to clarify the morphology 
of the spectrum and its dependence on two different quantum numbers. 



6.1 The Ground State 

Let us describe the properties of the ground state level with a major accuracy, starting from 
the result of its existence with a non- zero energy E , i.e. Eq = 19.831/? 2 . 



In FigJS] we plot the wave function \1/ 

|2 



in the 



u. 



0. 



plane and in FigJ6] the corresponding 

739.466. The 

knowledge of the ground state eigenfunction allows us to estimate the average values of the 



probability distribution whose normalization constant is equal to N 



e 2 = t 2 + 1 / 4 
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Table 2: The first ten energy eigen- Figure 4: xhe lines represen t i n a typical spec- 

values, ordered from the tral manner the first levels, 

lowest one. 

anisotropy variables u, v and the corresponding root mean square, i.e. 

(u) = -1/2 , (v) = 0.497 (46) 

Au = 0.266 , Av = 0.077 (47) 

This result tells us that in the ground state the Universe is not exactly an isotropic one, and it 
fluctuates around the line of the Misner plane /3_ = 0. However we observe that it remains lo- 
calized in the center of the Misner space far from the corner at v — > oo (the other two equivalent 
corners were cut out from our domain by the approximation we considered on the boundary 
conditions, but, because of the potential symmetry, they have to be unaccessible too). Thus 
we can conclude that the Universe, approaching the minimal energy configuration, conserves a 
certain degree of anisotropy but lives in the region where the latter can be treated as a small cor- 
rection to the full isotropy. Such a behavior is a consequence of the zero-point energy associated 
to the ground state which prevents the absence of oscillation modes concerning the anisotropy 
degrees of freedom; we numerically derived this feature but it can be inferred on the basis of 
general considerations about the Hamiltonian structure. In fact the Hamiltonian contains a 
term v 2 p 2 which has positive definite spectrum and cannot admit vanishing eigenvalue. 

6.2 Asymptotic Expansions 

In order to study the distribution of the highest energy levels, let us take into account the 
asymptotic behavior of the zeros for the modified Bessel functions of the third kind. 
We will discuss asymptotic regions of the plane (t, 0, n) in both the cases i) t 3> n and 
ii) t ~ n >> 1. 



i) For t>n the Bessel functions admit the following representation: 
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Figure 5: The ground state wave function Figure 6: We show the probability distribution 
is sketched i n correspondence to the ground state 

of the theory, according to the bound- 
ary conditions. 

where a = tc/A — y/t 2 — n 2 + t arccosh(t/n), p = n/t and Uk are the following polynomials 

ua(t) = 1 

«*-!(*) = - + 1 A 1 - M 2 )u k (t)dt (49) 



^ o JO 

Retaining in the above expression only those terms of o(j), the zeros are fixed by the following 
relation 

sin(^ - t + f(log(2) - log(p))) - ^ cos(^ - t + t(log(2) - log(p))) = (50) 



that in the limit for n/t <C 1, can be recast as follows 



In 



t\og(t/n) = Itt =4- t = (51) 

productlog( 



n 



where productlog(z) is a generalized function that gives the solution of the equation z = we w 
and for real and positive domain, it is a monotonic function of its argument. In (15"TI) / is an 
integer number that must be taken much greater than 1 in order to verify the initial assumptions 
n/t < 1. 

ii) In case the difference between 2n and t is o(n 1 / 3 )(t,n 3> 1), we can evaluate the first zeros 
by the expansion (worked out from formula (9) of [3] 



t = 2n + 0.030n 1/3 , (52) 

providing the lowest zero (and therefore the energy) for a fixed value of n and then the relation 
for the eigenvalues for high occupation numbers as 



2 

h 



An 2 + 0.12n 4/3 (53) 



In this region of the spectrum the condition t ~ 2n, reduces the 2 quantum numbers charac- 
terizing the system to a single one, the whole wave function becoming assigned by n. 

Ttip Fitinvp nnfllvsis slinws that as pffprt nf F)irirVilpt l-innnHarv rnnHitinn snrl in thp limit nf 



high occupation numbers, we get an analytical expressions for the discrete structure of the 
spectrum. In fact for large values of t it was possible to give analytical representations for the 
position of the zeros, but we emphasize that the request to deal with these approximations 
causes the loss of a large number of levels and prevents a complete discussion of the quantum 
chaos associated to the model. However the above expressions are of interest because allow 
to compare these results with the corresponding spectrum provided by Misner in his original 
work [TT] . Of course, such a comparison of the two results can take place only on a qualitative 
level; in fact, between Misner (a, (3 + , (3 J) and Misner-Chitre like (r, it, v) variables a crucial 
difference exists and it has to be recognized into the correspondingly different behavior of the 
potential walls. In the Misner scheme the domain available to the point Universe increases as 
a — > oo (a being the isotropic Misner variable) and, therefore, we deal with a non stationary 
infinite well; the Misner-Chitre like variables allow to fix the infinite potential walls into a 
time-independent configuration in the (it, 0, v) plane. However we can at least compare the 
behavior of the energy eigenvalues with respect to the occupation number n. 
In his original treatment, Misner replaces, for fixed a values (indeed he makes use of an adia- 
batic approximation, see [Ej), the triangular well by a square of equal measure and determines 
the energy spectrum in the form 

f.-y^TZ^K ,54) 

where S denotes the triangular well measure S = a 2 /A 2 , A being a numerical factor; above n± 

denote the occupation numbers relative to j3± respectively. 

In our approach, for sufficiently large n, we get the dominant behaviors 



E MC 



productlog(^) 



( 55 ) 

K [ii)2n + 0{n 2 / 3 ) 

Thus we see that, in case ii), a part from a numerical factor (which is different because of 
the different approximation made on the real domains), the only significant difference relies on 
the term a -1 ; in the case i) instead the situation is a bit different because 2 quantum numbers 
explicitly remain and we get a linear relation as far as we require I oc n by a factor much greater 
than the unity (indeed the function pro duct log (In/n) provides a smooth contributions in the 
considered region I 3> n). The difference of the factor a' 1 with respect to the Misner case can 
be easily accounted as soon as we observe that the following relations hold: 

(3± = ab±(u, v) , (56) 

where the functions b± can be calculated from the coordinates transformations 

1 + u + u 2 + v 2 
a = — e T — 

= -i + fr + a.'-^ (57) 

2y/3v 

(3 = -e^±^ 

' 2v 

but their form is not relevant in what follows. On the base of (|56|) . the measure of a domain D 
in the j3± plane reads 

/ dl3+d(3- = a 2 I \J(u,v)\dudv , (58) 

Jd JD' 

J being the Jacobian of the transformation associated to b±, while D' is the image of D onto 
the (it, 0, v) plane. As a consequence we see that between a measure s in the (3± plane and a 
similar one (even not exactly its image), there is a difference for a factor a 2 which immediately 

nrnvirlps an PTrnlanfltinn fnr thp HifFprpnpp in thp 9 pnprcrv sripptra 



7 The Inhomogeneous Picture and Conclusions 



At the end of our analysis, we wish to bring the reader's attention to some physical aspects of 
the inhomogeneous Mixmaster. 

First of all, the obtained dynamics regime is indeed a generic one; in fact, in a synchronous 
reference (for which dty a = 0), the inhomogeneous Mixmaster contains four independent (physi- 
cally) arbitrary functions of the spatial coordinates (y a (x l ) and E(x 1 )), available for the Cauchy 
data on a non-singular hyper-surface. 

Let us now come back to the full inhomogeneous problem, in order to understand the structure 
of the quantum space-time near the cosmological singularity. Since the spatial gradients of 
the configurational variables play no relevant dynamical role in the asymptotic limit r — > oo 
(indeed the spatial curvature simply behaves as a potential well) then the quantum evolution 
above outlined takes place independently in each space point and the total wave function of 
the Universe can be represented as follows 



where the product is (heuristically) taken over all the points of the spatial hypersurface. 
However, it is worth to recall that, in the spirit of the " long- wavelength approximation" adopted 
here, the physical meaning of a space point must be recovered on the notion of a cosmological 
horizon; in fact we are dealing with regions over which the inhomogeneity effects are negligible 
and this statement corresponds to super-horizon sized spatial gradients. On a classical point 
of view, this request is at the ground of the BKL approximation and it is well confirmed on 
the statistical level (see [4]); however on a quantum level it can acquire a precise meaning 
if we refer the dynamics to a kind of lattice space-time in which the spatial gradients of the 
configurational variables become real potential terms. In this respect it is important to observe 
that the geometry of the space-time is expected to acquire a discrete structure on the Planck 
scale and we believe that a regularization of our approach could arrive from a "loop quantum 
gravity" treatment [35] . 

Despite this local homogeneous framework of investigation, the appearance near the singular- 
ity of high spatial gradients and of a space-time foam (like outlined in the classical dynamics 
see [U |5]) can be easily recognized in the above quantum picture too. In fact the probability 
that in n different space points (horizons) the variables u and v take values within the same 
narrow interval, decrease with n as p n , p being the probability in a single point; in fact, these 
probabilities are all identical to each other and no interference phenomenon takes place. From 
a physical point of view, this simple consideration indicates that a smooth picture of the large 
scale Universe is forbidden on a probabilistic level and different causal regions are expected com- 
pletely disconnected from each other during their quantum evolution. Therefore, if we start 
with a strongly correlated initial wave function S (m, v) = S(r , u, v), its evolution toward 
the singularity induces increasingly irregular distributions, approaching (1591) in the asymptotic 
limit t — > oo. 

The main result of our presentation can be recognized in the clear correspondence established 
between the classical space-time foam and the quantum one. We have outlined how this link 
takes place naturally for a precise choice of the operator-ordering only and how the "energy 
spectrum" is a discrete one, due to the billiard structure of the point-like Hamiltonian. Fi- 
nally, we fixed as a new feature the zero-point "energy" for the ground state associated to the 
anisotropy degrees. 
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